Elucidating Fermi's Golden Rule via bound-to-bound transitions in a confined 

hydrogen atom 
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We demonstrate an effective method for calculating bound-to-continuum cross-sections by exam- 
ining transitions to bound states above the ionization energy that result from placing the system of 
interest within an infinite spherical well. Using photoionization of the hydrogen atom as an example, 
we demonstrate convergence between this approach for a large volume of confinement and an exact 
analytical alternate approach that uses energy-normalized continuum wavefunctions, which helps to 
elucidate the implementation of Fermi's Golden Rule. As the radius of confinement varies, the re- 
sulting changes in physical behavior of the system are presented and discussed. The photoionization 
cross-sections from a variety of atomic states with principal quantum number n are seen to obey 
particular scaling laws. 



I. INTRODUCTION 

While undergraduate students usually encounter 
bound-to-bound transitions,— they typically have little ex- 
perience with understanding and performing calculations 
on transitions of a bound-to- continuum nature, as exem- 
plified by ionization and molecular dissociation. Fermi's 
Golden Rule, 2 
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provides an expression for the transition rate from an ini- 
tial bound state i to the infinitude of continuum states 
{j} in a small energy range centered about W with p(W) 
being the density of states with the same symmetry as j. 
While this is often presented 3 for a time-independent per- 
turbation, H' , this rule can also apply for an sinusoidal 
perturbation, H'(t), with angular frequency w, by set- 
ting H' = H'(t)/(e iut + e - *"*), rendering H' again time- 
independent— ~& However, when applying Fermi's Golden 
Rule it is not immediately clear how to interpret the 
product of p(W), which is infinite for states in the con- 
tinuum where all energies are allowed, and the mean of 
the square of the matrix element which incorporates the 
space-normalized wavefunctions {0|}, which have zero 
amplitude as they extend over all space. 

In this paper we use two methods to resolve this math- 
ematical dilemma, and in doing so provide a concrete 
demonstration of the consistency between them. We 
thereby elucidate some of the subtleties of Fermi's Golden 
Rule which are often omitted in texts of introductory 
quantum mechanics jliLS but are important to under- 
stand for its applicability. The first method uses a par- 
ticular form of normalization for the continuum states 
called energy-normalization,— which we introduce below. 
The second method involves placing the system in an in- 
finite well (box normalization)^ which causes all states 
to be bound so that states in the previously continuous 
regime become discretized, allowing the familiar tech- 
niques from a bound-to-bound approach to be used to 
perform bound-to-continuum calculations. 



We use as an example photoionization of hydrogen, 
the process whereby the atom absorbs incident light such 
that the electron is no longer bound to the system. Per- 
forming calculations related to a confined hydrogen atom 
is not a new concept but past calculations of this 
kind have been done primarily to determine the proper- 
ties of the atom inside a confining structure, and here we 
use confinement as a tool to determine the properties of 
the unconfmed system. 

This paper is organized as follows: in Section |TT] we 
introduce the concept of a cross-section, apply it to more 
familiar bound-to-bound photoexcitations, and provide 
results specific to the hydrogen atom. In Section [TTT] we 
examine bound-to-continuum transitions, using an ana- 
lytical approach for the unconfmed atom, introduce en- 
ergy normalization, and reproduce, for later comparison, 
exact calculations of photoionization cross-sections for 
hydrogen, the archetypal charged binary system. In Sec- 
tion IIVI we implement box-normalization, demonstrate 
convergence for increasing radius of confinement, and 
show that calculations using this method agree with the 
exact calculations from the previous section. 

Two commonly encountered constants used through- 
out this paper are ao, the Bohr radius (= 4:Tre h 2 / (m e e 2 ) 
« 5.29xl0 _11 m), and Ry, the Rydberg unit of energy 
(=rn e e 4 /[(47re ) 2 2?i 2 ] « 2.18xl0~ 18 J). 



II. CROSS-SECTIONS FOR 
BOUND-TO-BOUND TRANSITIONS 

A physical parameter that quantifies the absorption of 
light by a sample of atoms with some distribution of ini- 
tial states is the cross-section, a. The average photon 
absorption rate per atom due to an incident narrowband 
photon flux, F (photons/unit time/unit area), is giveni^ 
by r = <j(lS) F, and thus the cross-section indicates the 
average effective area of an atom presented to the light at 
angular frequency w. By considering the photon absorp- 
tion rate of a single weak incident beam (i.e. assuming 
negligible population inversion) traveling in the positive 



2 



x-direction within a thickness dx of a sample, it can be 
shown that the intensity of the light, I(x), obeys 



dl , , 



(2) 
(2) 



where N is the number density of all particles. Eq 
is often given in textbooks as the definition for a. 

For bound-to-bound transitions between two partic- 
ular quantum states, a standard semiclassical treatment 
using the time-dependent Schrodinger equation for a nar- 
rowband, weak incident light beam gives* 9 
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where d is the electric dipole moment operator due to 
all charges in the atom, F e is the incident photon flux 
with polarization e which couples the transition and 
g(cj — ujq) is the lineshape of the transition with Jg du> = 
1. In general, however, many atomic levels are degenerate 
and light can couple many of these states, so the average 
absorption rate per atom for an ensemble of atoms is 
given by 



Ni/N 



(4) 



where the summation runs over all pairs of states coupled 
by the incident light, and Ni/N is the fraction of atoms 
in each initial state 

It is usual to report absorption cross-sections that as- 
sume the atoms are evenly distributed among all 2 J + 1 
possible M initial states. This distribution is frequently 
encountered, and T is then by symmetry independent of 
the polarization characteristics of the incident light. For 
ease in deriving this cross-section, we will arbitrarily con- 
sider incident light which is linearly polarized along the z- 
axis (7T-polarized, so AM=0). For a transition 7 J— y/J' 
this gives 



ttlj g(u — UJq) 
ce h 2J+1 



J2\(j'J'M\d z \jJM) 



M 



(5) 



where d z is the ^-component of the electric dipole mo- 
ment operator. For the particular case of nl^-n'i' tran- 
sitions within the hydrogen atom, neglecting electron and 
nuclear spin, this gives 

irue 2 g(u -uj ) i2 
0~ = ; — : — \Ur, 
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where the angular integral was evaluated using Ref. l20|, 
f mal is the greater of I and £', \£ — £'\ = 1, and the ra- 
dial matrix element is defined in terms of the well-known 
space- normalized radial wavefunctions, R n £, by 
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R* /e rRni r 2 dr. 
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It follows that the integrated cross-section, a commonly 
presented quantity that reflects the overall strength of a 
transition by integrating a over all incident frequencies, 
v, is given by 



o~ dv 
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III. CROSS-SECTIONS FOR 
BOUND-TO-CONTINUUM TRANSITIONS: 
ENERGY-NORMALIZATION 

One approach to evaluating the problematic product 
described in Section U is to use energy-normalized wave- 
functions for the continuum states, defined by>2i 



(9) 



Integrating over W on both sides for an arbitrarily small 
range centered about W gives another representation of 
this equation: 



W+AW 

J ^W'dW'dr = 1. 

W-AW 



(10) 



We examine a hydrogen atom with space-normalized 
continuum wavefunctions j = Y^Rwe (i.e. a ra- 
dial wavefunction with zero amplitude) and show that 
multiplication by \J p(W)' (which is infinite) gives 
energy-normalized finite-amplitude wavefunctions satis- 
fying Eq. (|10p . The left hand side of this equation takes 
the form 

W+AW 

VpYt*R*w?Jvp Y I n Rw>i dW dr 

W-AW 
W+AW 
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Rw'idndr — J r R we Ry/t dr = 1, 




where we have used in the first step that each p is dn/dW 
evaluated at W and W respectively, with W' bounded 
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by W ± AW, and in the final step the integral over n 
is unity since there is only one radial wavefunction Rw'i 
which is not orthogonal to Rwi- 

The photoionization cross-section for hydrogen can 
then be derived from Fermi's Golden Rule by including 
the density of states within the bra: (y/pj \ , and using for 
this energy-normalized continuum wavefunctions (energy 
W > 0) with radial component^ 
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The normalization constants herein can be derived by ap- 
plying Eq. (fT0|) to the asymptotic form of the radial wave- 
function, and these result in wavefunctions that when 
multiplied by r tends to an oscillatory function with con- 



stant amplitude of (aQ7r 2 H / Ry) 



-1/4 



for large 



Note 



that for W = — Ry /n 2 the functional form of Eq. (fl~2j) , 
which contains \F\ (a; 6; z), the confluent hypergeomet- 
ric function,— reduces to that of the well-known general- 
ized Laguerre polynomials for hydrogenic bound states. 

Again assuming a uniform distribution of initial m 
states and employing 7r-polarized light with electric field 
amplitude Eq, and identifying H' to be <I z Eq/2, Eq. ^ 
becomes 
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By equating Eq to 2u:F„h/ceo , the photoionization cross- 
section can be seen to be: 
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where fku = W + Ry/n 2 . By analogy to Equation ([7]), 
we define 



D nt ^ wt = I Rw£> r R n i r 2 dr, 



(15) 



which can be evaluated exactly for any n, £, and I' to 
determine the photoionization cross-section as a function 
of energy. We have calculated this cross-section from a 
variety of initial states and show them in the following 
section, demonstrating agreement with the result from 
box normalization calculations described below. 
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FIG. 1. (color online) The energies of states as a function of 
radius of confinement. The inset shows an additional curve 
representing the potential energy as a function of r on the 
same scale. 



IV. CROSS-SECTIONS FOR 
BOUND-TO-CONTINUUM TRANSITIONS: 
NORMALIZATION 



BOX 



A second approach commonly proposed to resolve the 
mathematical dilemma described in Section U is to place 
the system in a box of finite size, so that the problematic 
terms are neither zero nor infinite^ We now demonstrate 
this approach, showing that the results converge to the 
correct values as the volume tends to infinity. 

The potential energy of a hydrogen atom contained in 
a spherical well of radius with origin located at the 
atom's center of mass is 



U(r) = 



2R.y a 



r < r 
r > r . 



(16) 



We start by solving the Schrodinger equation to find 
the allowed energies W n ((r ) and energy eigenfunctions 
for such a confined hydrogen atom. As the system re- 
mains spherically symmetric, the wavefunctions retain 
the standard Yf 1 angular portion of the free hydrogen 
atom. 

We adopt the functional form of Eq. ([T2"|) for each 
confined radial wavefunction, d^w^e, and introduce the 
boundary condition: 



(17) 



Allowed energies under this condition were found nu- 
merically for fixed £ and tq using Mathematica, and 
agree with those found previously by Aquino et alJ^ The 
boundary condition represented by Eq. (|17l) has the effect 
of raising the energies from those of the free atom bound 
states, as shown in Figure Q] The inset illustrates that 
the energies begin to rise more dramatically when the ra- 
dius of confinement, tq, is such that the state no longer 
has a region of exponential decay. States previously in 
the continuum no longer exist since there is no contin- 
uum; we call those bound states which rise to an energy 
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greater than the ionization limit pseudocontinuum states, 
and relabel their energies W, dropping the subscript n£ 
to distinguish them from those of the (Coulombically) 
bound states. 

These bound and pseudocontinuum wavefunctions 
were then numerically space-normalized in the usual way. 
We wish to find prefactors for each of the initial and fi- 
nal states which will have the property of causing the the 
radial matrix element to tend, as rg— >oo, to that given 
in Section UTT1 for the free atom. 

We define the density of final states with a par- 
ticular symmetry p (W) for the confined atom to be 
2/(W + —W~) where W + (W~) is the energy of the adja- 
cent state with the same symmetry above (below) |9tvt^)- 
Thus, as ro goes to infinity so too will the density of states 
above the ionization limit. Furthermore, from Section Hill 
as ro — > oo, pseudocontinuum states 9iw£ ( r ) assume zero 
amplitude in such a way that 



y/p(W)'<Rwe(r) ^ R% t (r) 



(18) 



and for (Coulombically) bound states, trivially 

Kw ne e{r) -> R n e(r) , (19) 



where in each case CH, the confined-atom wavefunc- 
tion, is space- normalized, and R^e an d Rnt are energy- 
normalized and space-normalized, respectively. The 
free atom bound-to-continuum radial matrix element in 
Eq. (|15j) can therefore be readily found from 



I Vp(W)'K* wl , m Wnl t r 2 dr (20) 
-> D n £^wi'- (21) 
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We have calculated the photoionization cross-sections 
for the hydrogen atom using the infinite well approach 
for Is —¥ Wp for a variety of values of ro and show in 
Fig. [5] that it converges towards those found using ex- 
act calculations from Section HlH which agree with those 
from Ref. [H. The data at r = 20a and r = 50a 
is sparse as the densities of states at these radii of con- 
finement are quite low; all allowed final energies in the 
range shown are plotted. With r = 1500ao, there 
are sufficiently many states in the energy range shown 
to produce a smooth curve on the scale of the figure. 
The consistency with Fermi's Golden Rule is now ap- 
parent, in which the density of states evolves from an 
inverse energy difference: multiplication of the square of 
the matrix element by the density of states is mathe- 
matically equivalent to instead using an energy- rather 
than space-normalized continuum wavefunction. It can 
be seen that the cross-section converges to the true con- 
tinuum case from above as ro increases. This can be 
attributed to the fact that while r'ZKwi tends to a con- 
stant amplitude as r —¥ oo, it does so from below; thus at 
low r the amplitude is slightly reduced from this value. 
Therefore when the wavefunction is truncated at finite ro 
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FIG. 2. (color online) The photoionization cross-section for 
Is — > Wp. The results from the infinite well approach are 
shown with triangles for ro = 20ao , circles for ro = 50ao , and 
a dashed line for ro = 1500ao, and the solid line is the ex- 
act calculation. The last two are shown in the inset for large 
energies, along with the dotted line showing the Born approx- 
imation. 



before reaching its final amplitude, 9ty^ decreases dis- 
proportionately to increases in ro. This results in the 
Dw„ e ^-we' of Eq. ([20]) taking slightly too large a value for 
finite ro. The inset of Fig. [2] shows that this calculation 
works for a large range of final-state energies, over the 
entire non-relativistic regime. At high energies, the cross 
section can be seen to approach that predicted by the 
Born approximation, 25 which for a transition ns — > Wp 
takes the form: 



W- 
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(22) 



An alternative measure of the likelihood of a 
bound-to-continuum transition is the oscillator strength 
distribution^ which is related to the cross-section, 
by 



ch 



df_ 

dW 8agRyV 



(23) 



and for completeness these values are given in the figures 
as well. 

We also show in Fig. [3] a plot to illustrate that this ap- 
proach works for a range of values of initial n. Higher n 
results in a greater cross-section for low ejection energies, 
and a lesser cross-section for high ejection energies. This 
relates to the general observation 27 that high n (Ryd- 
berg) states are strong absorbers of long-wavelength ra- 
diation. The peak cross section occurs at threshold; for 
ns — > (W=0)p, this simplifies to 



2 Wa 2 fj^ 2- (l-n)< m > L^j_ m {-2n) 

\m=0 
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(24) 

where iffi (x) is a generalized Laguerre polynomial, and 
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(1 — n) (2 — n) • • • (m — n) is the Pochham- 
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FIG. 3. (color online) The photoionization cross-section for 
{5, 10, 20} s —> Wp. The exact calculations are shown as lines, 
with symbols showing results using the infinite well approach 
with ro = 1500ao. Ry is the Rydberg unit of energy. 



mer function, which becomes zero for any m > n, forcing 
the sum to terminate. 



V. CONCLUSION 

Wc have presented two methods for calculating bound- 
to-continuum cross-sections: the first employing energy- 
normalization and the second, box-normalization, and 
demonstrated convergence between them, using pho- 
toionization of the hydrogen atom as a concrete exam- 
ple. We thereby elucidate some quantitative aspects of 
Fermi's Golden Rule and show the product of the density 
of states and a space-normalized continuum wavefunction 
to be mathematically equivalent to using only an energy- 
normalized continuum wavefunction. 
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